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Abstract 

We analyze the semileptonic Dq Kiiu transition with q = u, d, s, in the 
framework of the three-point QCD sum rules and the nonleptonic D — > Kiir decay 
within the QCD factorization approach. We study Dq to i<'i(1270) and iiri(1400) 
transition form factors by separating the mixture of the -ftri(1270) and i^i(1400) 
states. Using the transition form factors of the D Ki, we analyze the nonleptonic 
D — > i^ivr decay. We also present the decay amplitude and decay width of these 
decays in terms of the transition form factors. The branching ratios of these channel 
modes are also calculated at different values of the mixing angle 6ki and compared 
with the existing experimental data for the nonleptonic case. 
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1 Introduction 



Analyzing the semileptonic decays of the charmed Dg mesons is very useful for determina- 
tion of the elements of the Cabibbo-Kabayashi-Maskawa (CKM) matrix and also leptonic 
decay constants of the initial and final meson states. The semileptonic Dg — > Kiii/ tran- 
sition could give useful information about the internal structure of the Dg meson. Investi- 
gating the nonleptonic decays such as D ^ Kin can also be important for interpretation 
of the structure of the lightest scaler mesons [1] . 

Prom the experimental view, the physical states iri(1270) and ^1(1400) are the mixtures 
of the strange members of two axial-vector SU{3) octets l^Pi(i^i^) and 1^Pi{Kib). The 
and Kib are not mass eigenstates and they can be mixed together due to the nonstrange 
light quark mass difference. Their relations with the i^i(1270) and i^i(1400) states can be 
written as [2-4]: 

I ii'i(1270) > = I KiA > sin0K, + | Kib > cosOk^, 

I KiilAm) > = \KiA> cosOk^ -\KiB> sinOK,. (1) 

The angle has been obtained with two-fold ambiguity | |~ 33°, as given in Ref 
[3]. Also in Ref [6] 35° <| 6ki |< 55° has been found. In this paper we use Ok^ in the 
interval 37° <| |< 58° [4, 7]. The sign ambiguity for 6ki is due to the fact that one can 
add arbitrary phases to | Kia > and | Kib > states. 

The QCD sum rules approach has been successfully apphed to a wide variety of problems 
in charm meson decays. The semileptonic decays Ds — > foii^, Dg — > [1], D ^ K iu [8], 
D+ K^*e+Ue [9], D Tiiu [10], D piu [11], Df [12] and D K*iu [13] have 

been studied in the framework of the three-point QCD sum rules. As a nonperturbative 
method, the QCD sum rules has been of interest and it is a well established technique in 
the hadron physics since it is based on the fundamental QCD Lagrangian (for details about 
the QCD sum rules approach see for instance [14]). 

In the present work, wc study the semileptonic decays of the Dg Kiiu in the frame- 
work of the three-point QCD sum rules. The long distance dynamics of such transitions 
can be parameterized in terms of some form factors calculating of which play fundamental 
role in the analyzing of such type transitions. Considering the contributions of the oper- 
ators with mass dimension d = 3, 4, 5 as condensate and non-perturbative contributions, 
first we calculate the transition form factors of the semileptonic Dg Kiii'[q = u, d, s) 
decays. Using these form factors, the total decay width as well as the branching ratio for 
the aforementioned transitions are also evaluated at different values of the mixing angle. 
Having computed the form factors of the D ^ Ki, the amplitude and decay rate of the 
nonleptonic Du,d Kitt decays are also computed in terms of those form factors using the 
QCD factorization method (for more about the method see [15-17] and references therein). 

The paper is organized as follows. The calculation of the sum rules for the relevant 
form factors arc presented in section2. In calculating the form factors, first wc consider the 
general {Ki\ state. Then, using the definition of the G-parity conserving decay constant 
< I Jria I ^ia{p':£) >— Jkia'^Kia^^ G-parity violating decay constant < | J^^^ | 
KiB{p',e) >— fK^gj_{i GeV)aQ^^'^mj(Be'', where a||'^^^ is the zeroth Gegenbauer moment 
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of KiB state and it is zero in the SU{3) symmetry limit, we obtain the form factors of the 
D KiAi^B) states. Finahy, considering Eq. (1), we separate the (f^i [1270(1400)] | states 
and derive form factors of the D — > i^i [1270(1400)] transitions. The decay rate formulas 
for semileptonic and nonleptonic cases are presented in sectionS. We derive the decay rate 
formula for D — > Xitt decay using the QCD factorization method in tree level. Section 4 
is devoted to the numeric analysis of the form factors as well as the branching fractions of 
the considered semileptonic and non-leptonic decays at different values of the mixing angle, 
and discussions. A comparison of our results for the branching ratios for the non-leptonic 
case with the existing experimental data is also made in this section. 

2 Sum rules for Dq Kiiu transition form factors 

The Dq Kiiu with q — u,d,s decay governed by the tree level c ^ q' {q' — d, s) 
transition (see Fig .1). 




u d s 



(1) (2) (3) 



Figure 1: Semileptonic decays of Dq to Ki. Diagrams 1, 2 and 3 are related to the D^ — > 
K{tv, D+ K^eiy and Df K^iu, respectively. 

In the standard model, the effective Hamiltonian responsible for these transitions is 
given as: 

G 1 / 
^eff = -^c<{ ^ 7m(1 - 75) I i 7m(1 - 75)c, (2) 

where, is the Fermi constant and are the CKM matrix elements. The decay ampli- 
tude for Dq — > K\iv is obtained by inserting Eq. (2) between the initial and final meson 
states. 

M = ^K,' ^ 7m(1 - 75)/ < i^i(p', e) I i 7^,(1 - 75)c | Dq{v) > ■ (3) 

The next step is to calculate the matrix element appearing in Eq. (3). Both axial 
and vector parts of the transition current give contribution to this matrix element and it 
can be parametrized in terms of some form factors using the Lorentz invariance and parity 
conservation as follows: 

< K^ip', e) I ^'7,.75C I Dq{p) - '^^'U ^.a/?£>V^ (4) 



(5) 



In order for the calculations to be simple, the following redefinitions are used 



^ V? J - 7- — , ^ N > ^0 [Q ) - h 



(mo, + ) (mc, + ) 

where the Fy'^^^q^), Fo'^^'iq^), F^'^^^q^) and F^'^^'iq'^) are the new transition 
form factors, P/^ — {p+p')ix, q^ — ip—p')ij. and e is the four-polarization vector of the axial 
Ki meson. 

Based on the general philosophy of the three-point QCD sum rules technique, the above 
form factors in Eq. (6) can be evaluated from the time ordered product of the following 
three currents. 

^%-^\p',p", e) = ^'J d'xd'ye+^^'^-^^y (o |t { Ji.,.(a:)jf (0)4^(1/) }| o) , (7) 

where, Jxivix) — qi'Jv'JbS {qi — u,d) , Joqiu) — TibC are the interpolating currents of the 
and Dq and = qinC and — q^ii^^c are the vector and axial-vector parts of the 
transition current, respectively. 

The above correlation function is calculated in two different approaches: On the quark 
level, it describes a meson as quarks and gluons interacting in a QCD vacuum. This is called 
the theoretical or QCD side. In the phenomenological or physical side, it is saturated by a 
tower of mesons with the same quantum numbers as the interpolating currents. The form 
factors are determined by matching these two different representations of the correlation 
function and applying double Borel transformation with respect to the momentum of the 
initial and final meson states to suppress the contribution coming from the higher states 
and continuum. We can express the correlation function in both sides in terms of four 
independent Lorentz structures: 

n£^"''^ = e^vo^fi py^Iiv + 9^.uI^o + P^.PuI^l + q^.PuI^2 ■ (8) 

To find the sum rules for the related form factors, we will match the coefficients of the 
corresponding structures from both representations of the correlation function. 

First, we calculate the aforementioned correlation function in the phenomenological 
representation. Inserting two complete sets of intermediate states with the same quantum 
number as the currents Jk^ and Jo^ to Eq. (7), we obtain 

< I Jk,. I Mp\e) X K,{p',e) \ Jj^-^ \ Dgjp) >< Dqjp) | 4 J > ^ 

the higher resonances and continuum. (9) 
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In Eq. (9), the vacuum to initial and final meson states matrix elements are defined as: 



.2 



f 2 

< I J^J K,ip') >= fK.mK.e" , <0\Jd,\ D,{p) z^^^, (10) 

TTIq n~ ^^^^ 

where fxi and are the leptonic decay constants of Ki and Dg mesons, respectively. 
Using Eq. (4), Eq. (5) and Eq. (10) in Eq. (9) and performing summation over the 
polarization vector of the Ki meson, we get the following result for the physical part: 

n,. (p,P,g) - -(^^ + ^^) (p'2_„,2^j(p2_^2^^) x[^o {<l)9,. 

+ excited states. (11) 

The coefficients of the Lorentz structures ie,j,i,ai3P°'p^ , Qixv, P^Pu and Qf^p^ in the correlation 
function H^""^ will be chosen in determination of the form factors F^'-''^^^^ {(f), Fq'^"-'^^^ (g^) , 
^^^(.)-^i(g2) F^^'^^'^'ig^), respectively. 

On the QCD or theoretical side, the correlation function is calculated in the quark and 
gluon languages by the help of the operator product expansion (OPE) in the deep Euclidean 
region where p^ -C {rric + r«g)^, p -C (m^ + m\). In Eq. (7), using the expansion of the 



time ordered products of the currents, the three-point correlation function is written in 
terms of the series of local operators with increasing dimension as the following form [18]: 

+ (C5)M.*<^a/3G"^* + (C6)^,^r W'* , (12) 

where. Gap is the gluon field strength tensor, {Ci)^^, are the Wilson coefficients, / is the 
unit matrix, is the local field operator of the light quarks, and F and F are the matrices 
appearing in the calculations. Taking into account the vacuum expectation value of the 
OPE, the expansion of the correlation function in terms of the local operators is written as 
follows: 

^i^ApIpI q') = Com. + C3,u(^^) + C,^,{G^) + C^^^i^a^^G'^P^) 

+ C6^,(^FW'*) . (13) 

In Eq.(13), the contributions of the perturbative and condensate terms of dimension 
3,4, and 5 as non-perturbative parts are considered. The diagrams for the contributions of 
the non-perturbative part are depicted in Figs. 2, 3 and 4. It's found that the heavy quark 
condensate contributions are suppressed by inverse of the heavy quark mass and can be 
safely removed (see diagrams 4, 5, 6 in Fig. 2). The light q quark condensate contributions 
are zero after applying the double Borel transformation with respect to both variables p^ 
and p since only one variable appears in the denominator (see diagrams 1, 2, 3 in Fig. 2). 

Our calculations show that in this case, the two-gluon condensate contributions (see 
diagrams in Fig. 3) are very small in comparison with the quark condensate contributions 
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Figure 2: The quark condensate diagrams without any gluon and with one gluon emission, 
and we can easily ignore their contributions in our calculations. 



Therefore, the main contribution in the non-perturbative part comes from the q-quark 
condensates, (see Fig. 4). 

As a result, in the lowest order of the perturbation theory, the three-point correlation 
function receives a contribution from the perturbative part ( bare-loop contributions of 
diagrams in Fig. 1) and nonp-erturbative part (contributions of the diagrams shown in Fig. 
4) i.e., 



Mp',p \ q') = nr (p^p", q') + nr-^-(/,p", q') . (14) 

Using the double dispersion representation, the bare-loop contribution is determined: 

1 r r i:^"'{s s' a^) 
nf = -- — - / / , ' ' \' \^ dsds' + subtraction terms , (15) 
(27r)^ J J [s — p^) [s' — p'^) 

The following inequality is responsible for obtaining the integration limits in Eq. (15). 

2ss' -\- {s-\- s' — q^){ml — ml — s) + 2s{ml — m^,) 
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Figure 3: Diagrams for two-gluon condensate contributions. 



Yu(1-Y5) 




(1) 




(2) 




(3) 



Figure 4: Diagrams for q-quark condensates contributions. 



where A (a, 6, c) = + 6^ + — 2ah — 2ac — 2bc is the usual triangle function. 

By the help of the Cutkosky rule, i.e., replacing the propagators with the Dirac-delta 
functions: 



^2 _ 

the spectral densities p^^^{s, s', q^) are found as: 

pv = 4:Nc Io{s, s', q^) {Bi{mc - niq) - .62(^5/ + ruq) - mq} , 



(17) 



Po = -2A^c s',g ) {A(mq + rry ) - A'(mc - mg) - 4/li(mc - mg) 
+2im?{mc — rriq — m ') + mq{2mcm > — u)} , 



pi = 2Nc Io{s,s',q ){Bi{mc-3mq) - B2{mq + mg') + 2A2{mc-'mq) 
+2A3{mc - rUq) - rUq} , 
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+B2{mq + 171^')+ mq} . 



where 



(18) 



Io{s,s',q ) — , ,„ . — ^, 

A(s,s',g2) = + s'^ + - 2sg^ - 2sV - 2ss', 

A2 = — -^^ — ^[2A'V + 6AV'-8m2ss''-6AAVxi 

+A'^M^ + 2m^sV], 
A3 = — — -[-3AW-3A'\s + Wms's + 4AA'ss' 
+2AAV-my], 



where, m = s + s' — g^, A = s + — , A' = s' + — m^, and iVc = 3 is the color 
factor. 

The corresponding non-perturbative part of the considered structures are obtained as 
follows: 

rrnon-per/ 2 '2 2n ^ fl rnqUlq' 1 777,^777.^ Ulqi'^mq^ 1 TJlqi'^mQ^ 

I 2 rr 2 r r rr 2 rr 

1 rr, 2^ 2 1 ^„2^ 2 1 ^ ,2^ 2 1 ^ ,2^„2 i ^2^ 2 



2 




3 r^r'^ 




2 rV^ 


3 r^r' 


1 


2 2 

mo q 


2 2 

rric rriq 


1 


2 2 1 

rriQ rric 1 


niQ^mcmqi 


3 






2" 







2 + 2 ^2^/2 
1 mn^ 1 

+ :t 



3 r^r' 



(19) 



/j^ /y/y^ ^ ry'2iryf ^ iytyf ^ iy'2,iyf 



I /I > ' y ' ' /I /yZj/yi /] ly ly • 



1 22i 22i 24i 24i 24 

1 niQ niqi 1 mo q 1 mo mc 1 mo mq' 1 mo q 
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1 24i 42i 22i 22i22 

1 TJiQ rriqi I rric rriq I rric rriq I rriq/ rriq 1 q rriq 

1 rriqi'^m^ 1 rUqi'^m^ 1 mo^ 3 mo^mcmqi ^ 1 mo^m^mqi'^ 

1 mQ^nii^q^ 1 niQ^mqi'^if' 1 mqiTif^mqi 1 mqiricmqi'^ 

1 rnqTricq^ 1 rnqniqiq^ 1 mqimc^mqi 1 rnqTricmq''^ 1 rrig^ 
^ j^2j^/ ^ rr'^ 4 rr'^ 2 rr'^ 2 rr' 

1 m^m^ 1 q^m^ ^ 1 (^mq 1 mcvaq ^ 1 mcva^ 

1 mcmqimq^ ^ 1 mcmqimq^ 1 mcmqimo^ ^ m^tTiqimq^ 

1 mc^mqimo^ 1 rric^mqi'^mq^ 1 m^mqi'^mQ^ 1 rrii^q^mq^ 

1 rriQ^'m^if' 1 m^q^m^ 1 mcmqi(fmq^ 1 mcmgiq^mo^ 

1 niqi'^q^m^ 1 m^mqi'^mq^ ^ 1 m^mqi'^mQ^ mcmqi^mq^ 

1 mcmqi^mo^ 1 mqi^q^mq^ 1 mo^mqi'^q^ 1 mc^mq'niq'^ 
~2 rr^^ ^ 2^ rr^^ 4 rr^^ ^ 2^ rV^ 
1 mc^mg/mo^ 1 m^mqi'^m^ 1 mcmqi^m^ 1 mcmg/^mo^ 1 



(20) 



rrnon-per/ 2 '2 2n ^ _ ^ / 1 mgmg' 1 m^mc , 1 Ulq^'^mq^ 1 Ulqi'^mo^ 

Hi Mp,P ,0 = I - 4^:^:7^ + 47^+2 ^^/3 -4 ^^;3 

1 rric^mq^ 1 mo^m^ 1 niqi'^mq^ 1 mg/^mo^ 1 (fniq^ 

1 rriQ^q^ 1 m^rrif^ 1 niQ^m^ 1 m^^ 1 

1 mo^mcmqi \ 



2^2 



12 r2r 



nnon-per/ 2 '2 2n _ . - . j 1 ^g^g' 1 ^g^c 1 rUg^'^mq^ 1 mg/^mp^ 



(21) 



/j^ lyiyf ^ y'2iy/ 2 lyiyf'^ ^ lyiyf' 

1 m^mg^ 1 mo^m^ 1 m^/^mg^ 1 mqi'^mo'^ 1 g^mg 
+ -_ 



1 mo^?^ 1 rricmg^ 1 rriQ^m^ 1 m^^ 1 mo^ 



^ y'2:yf'^ 2 y^3i^/ ^ ry^yf 2 y^yf 2 y'^yf 

1 mo^mcn 
12 rV^ 



imj^l (22) 
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where r = — , r' — p' — m^, . 

Equating two representations of the correlation function and applying the double Borel 
transformation using 

™2 



1 ,m i-^r e ^ 



p^-ml' r(m) (Mf)"^' 
,2^, 1 (-1)" 



the sum rules for the form factors F^^^^~'^^ are obtained as: 

-Fi = --^ n „ e ^1 e - — — / ds / s , ? je^^i e'^a 



+ Ml Ml B^Ml) B^,,{Ml) [nr-^^^(/,p",?^)]|, (24) 



where i = V^, 0, 1 and 2, Sq and ,Sq are the continuum thresholds in pseudoscalar Dg and 
axial-vector Ki channels, respectively and the lower limit in the integration over s is as 
follows: 

_ (m^ + q'^-ml- s'){mls' - my) 
[m^ — q )[m^ — s') 

In Eq. (24), to subtract the contributions of the higher states and the continuum the quark- 
hadron duality assumption is also used, i.e., it is assumed that 

phigherstates^^^ s') = p°^^(s, s')e{s - So)e{s - s'^). (26) 

Here, we should stress that in the three-point sum rules with double dispersion relation, 
the subtraction of the continuum states and the quark-hadron duality is highly nontrivial. 

For > values, their may be an inconsistency between double dispersion integrals in 
Eq. (24) and corresponding coefficients of the structures in the Feynman amplitudes in 
the bare-loop diagram. In this case, the double spectral density receives contributions 
beyond the contributions coming from the Landau- type singularities. This problem has 
been widely discussed in [9]. Here, we neglect such contributions since with the above 
continuum subtraction and the selecting integration region the contribution of the non- 
Landau singularities is very small comparing the Landau type singularity contributions. 

Now, as wc mentioned in the introduction section, the '' ^^^^^ form factors arc ob- 
tained from the above equation by replacing fx^ by the G-parity conserving decay constant 
fxiA ^^'i G-parity violating decay constant fx^g = fK^^j_{^ GeV)aQ^^'^ and mK^ with 

1A(S) ' 



IS 



Fi " ' e'^ie J- ds / dspi{s,s ,q^)e'^ie'^2 

+ Ml Ml BAMl) B^,,{Ml) [Hr-^^^(/,p",g^)]|. (27) 
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Also, using Eqs. (1, 4, 5, 6), the form factors of the y.^9^-'^i[i2™(i400)] found as follows: 

j,Dg^Ki(1270) 
J1,2,V 

Jo 





+ rriKi 






mo. 


+ friK.A 


mo. 




mo. 




, mo. 







1 /v-w-i ' 


mo. 
















, mo, 





(28) 



3 Decay amplitudes and decay widths 
semileptonic 

Using the amplitude in Eq. (3) and definitions of the form factors, the differential decay 
widths for the process Dq — > Kilv are found as follows: 

2 



cq 



1927r3 



D„ 



cq 



dq^ 



1927r3m|) 



q'X'/'{ml,ml^,q')\H,f , 



(29) 



where, 

H±{q') = 

Ho{q') = 



A ' [jU'n ■ nil- ■ <r) 

{mn, + mK,)Mq') T \ ^ Mq') , 

rriDg + rriKi 



2mKiVq^ 



[m 



2 2 2x/ , ^./2^ K^h.^rnl^A 



D, -niK,-q ){mD, + mK^)fo{q ) 



The ±, in the above relations belong to the Ki helicities. The total differential decay 
width can be written as 



dq^ 



dq^ 



+ 



dq^ 



(30) 



where. 



dq^ 
dq^ 



dq'^ ' 
dr+{Dg^Kiiu) dT_{Dg^ Kiiu 



dq"" 



+ 



dq"" 



(31) 



and (d^) is the longitudinal (transverse) component of the differential decay width. 
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nonleptonic 



In this part, we study the decay amphtude and decay width for the nonleptonic D — > Xitt 
decay. The effective Hamiltonian for this decay at the quark level is given by (see for 
example [19] and references therein): 

Heff = ^{V^sVUC.O.-rCO,)]. (32) 

Here Oi and O2 are quark operators and they are given as: 

Oi = {SiCi)v-A{Ujdj)v-A, O2 = {SiCj)v-A{Ujdi)v-A: (33) 

where {qiq2)v±A = gi7^(l ± 75)52- 

The Wilson coefficients Ci and C2 have been calculated in different schemes [20]. In the 
present work, we will use Ci(mc) = 1.263 and C2(mc) = —0.513 obtained at the leading 
order in renormalization group improved perturbation theory at // ~ 1.3 GeV [21]. 

Now, we calculate the amplitude Aior D Kitt decay. Using the factorization method 
and definition of the related matrix elements in terms of the form factors fy^^^ , fo^^^ , fi^^^ 
and /^~*^^ in Eqs. (4-6), we obtain this amplitude as follows: 

A''-''^- = ^ «i} A [i^^^^^'^K)], (34) 

where, 

F^-^-(m^) = [{nin + mK,)fo{rnl)-{mn-mK,)h{ml)- , ^'^'^ , ^^]-(35) 

yrriD + rriKi) 

The e stands for polarization of Ki, p is four momentum of D, is the pion decay constant, 
ai = Ci + -^C2 and is the number of colors in QCD. 

Now, we can calculate the decay width for D Kin decay. The explicit expression for 
decay width is given as follow: 

T(D ^ Kin) = , \Vj^\Vud\^ al fl 

Ximl,ml^,ml)l [F^-^-(m^)]^. (36) 



4 Numerical analysis 

From the sum rules expressions of the form factors, it is clear that the main input parameters 
entering the expressions are condensates, elements of the CKM matrix V^^/ , leptonic decay 
constants foq, Jk^a and fK^g±: Borel parameters and M| as well as the continuum 
thresholds Sq and Sq. We choose the values of the condensates (at a fixed renormalization 
scale of about 1 GeV) , leptonic decay constants , CKM matrix elements, quark and meson 
masses as: < uu >=< dd >= -(0.240 ± 0.010 GeV)^, < ss >= (0.8 ± 0.2) < uu >, 
ml = 0.8 ±0.2 GeV^ [22], | Vcs \= 0.957 ± 0.110, | |= 0.230 ± 0.011 [23], /do = 
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fD± = 0.222 ± 0.016 GeV [24], /d, = 0.274 ± 0.013 GeV [25], fx,^ = 0.250 ± 0.013 GeV, 
/^^^^ = 0.190±0.010 GeV [2], m„(l GeV) = (1.5-3.3) MeV, md{l GeV) = (3.5-6) MeV, 
ms{l GeV) = (104:1:11) MeV, = 1.27;^! GeV, moo = 1.864 GeV, mD± = 1-869 GeV, 
mo, = 1.968 GeV, mi^,(1270) = 1.27 GeV, mi^,(1400) = 1.40 ^€^[23], rriK^^ = 1.31 ± 
0.06 GeV, niK.s = 1-34 ± 0.08 GeV and al'^'" = -0.19 ± 0.07 [2]. 

The sum rules for the form factors contain also four auxiliary parameters: Borel mass 
squares Mf and M| and continuum thresholds Sq and Sq. These arc not physical quantities, 
so the form factors as physical quantities should be independent of them. The parameters 
So and s'q, which are the continuum thresholds of Dg and Ki mesons, respectively, are 
determined from the condition that guarantees the sum rules to practically be stable in the 
allowed regions for and M|. The values of the continuum thresholds calculated from the 
two-point QCD sum rules are taken to be sq = (6 — 8) GeV^ and s'q = (4 — 6) GeV'^. The 
working regions for Mf and M| are determined requiring that not only the contributions 
of the higher states and continuum are small, but the contributions of the operators with 
higher dimensions are also small. Both conditions are satisfied in the regions 4 GeV"^ < 
Ml < 10 GeV^ and 3 GeV^ < M| < 8 GeV^. 

The values of the form factors at = are shown in Tables 1 and 2. Note that, the 
values of the /j(0) for D'^ — ^ K^iu and D"^ —>■ K^iu are approximately equal, so the values 
in Table. 1 refer to both decays. 





37 58 -37 -58 




37 58 -37 -58 


nD^Ki{1270) 

Jv 


3.19 1.82 4.00 2.95 


Jv 


-3.37 -4.34 2.27 3.60 


Jo 


-0.74 -0.42 -0.93 -0.68 


Jo 


0.72 0.92 -0.49 -0.77 


rD^Ki{1270) 
Jl 


0.34 0.19 0.44 0.34 


rD^Ki{1400) 
Jl 


-0.38 -0.49 0.23 0.38 


nD^Ki{1270) 
J2 


2.56 1.46 3.24 2.36 


^D^_ft:i(1400) 
J2 


-2.70 -3.49 1.82 2.90 



Table 1: The = values of the form factors of the D Ki£u decay for = 8 GeV^, 
M| = 6 GeV"^ at different values of 9ki ■ 

The dependence of the //^''^^^^(O) on Ok-^ at = is depicted in Figs. 5-8, in the 
interval —58° < 9ki < 58°. In Figs. 6 and 8, as it is seen, all of the form factors con- 
tact at one point. Also each form factor in Figs. 5 and 7, has one extremum point. 
These extrema as well as the contact points have been specified in Figs. 5-8. It is in- 
teresting that in the Dg A'i(1270)£z/ and Dg Ki{lAOO)ii' cases, the extrema and 
contact points of the form factors are nearly at —8°. The sum rules for the form factors 
are truncated at about = 0.15 GeV^ and = 0.25 GeV^ for q = u{d) and s cases 
of the Xi(1270), respectively. These points for Xi(1400) state are q'^ = 0.22 GeV^ and 
q'^ — 0.32 GeV'^ for u(d) and s cases, respectively. To extend the results to the full physical 
region, i.e., < < (^d^ "~ i^Ki)'^ GeV"^, we look for a parametrization such that: 1) this 
parametrization coincides well with the sum rules predictions below the points at which the 
form factors are truncated and 2) the parametrization provides an extrapolation to q^ > the 
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37 58 -37 -58 




37 58 -37 -58 


Jv 


3.90 2.22 4.86 3.58 


j.D+^K0{1400} 

Jv 


-4.09 -5.27 2.76 4.40 


Jo 


-1.15 -0.65 -1.44 -1.07 


„D+^iC<'(1400) 

Jo 


1.12 1.44 -0.76 -1.20 


j,D+-*KO (1270) 
Jl 


-0.54 -0.31 -0.66 -0.50 


Jl 


0.57 0.73 -0.39 -0.61 


rD+^K0{1270) 
J2 


5.89 3.36 7.33 5.40 


„D+^A-0(1400) 
J2 


-6.19 -7.97 4.18 6.64 



Table 2: The = values of the form factors of the Ds Ki£u decay for Mf = 8 GeV^, 
M| = 6 GeV^ at different values of 9ki ■ 

truncated points, which is consistent with the expected analytical properties of the form 
factors and reproduces the lowest-lying resonance (pole). This resonance in the Dg channel 
is D*{J^ = 1^) state. Following references [26, 27], which describe this point in details, we 
choose the following theoretically more reliable fit parametrization: 

1 — — 2 — J- — — 2 — 

The values of the parameters a, h and mfn are given in Tables 3-6 at different values of 
the mixing angle O^^- From this parametrization, we see that the itld* pole exist outside 
the allowed physical region and related to that, one can calculate the hadronic parameters 
such as the coupling constant Qdd'Ki (see [28,29]). 
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a b rrifit 




a b mfit 


fv (r) 


o oo A OA 1 or 

6.Q6 -U.d4 1.25 


fv \t) 


C f\A O CTH' 1 OCT 

-5.94 2.57 1.25 


.fo ir) 


-2.05 1.31 1.3(:) 


fo {<]) 


2.04 -1.32 1.36 


fi (r) 


0.4o -0.12 1.27 


fi (r) 


n m A oi 1 0*7 

-0.59 0.21 1.27 


/2 (r) 


2.97 -0.41 1.29 


f2 (r) 


o 1 /I o /I /I 1 oo 

-3.14 0.44 1.29 




4.08 -0.18 1.28 


D+ -*K°{1A00) / 2\ 

fv ir) 


'T O'T O '^O 1 OO 

-7.87 3.78 1.28 


j,D+ ^K0{1270),2\ 
Jo [Q ) 


-o.oD Z.41 l.oi 


j,D+^KO{UOO), 2\ 

Jo [Q ) 


o.OD -1.94 l.ol 


j,D+^KO (1270) ^^2^ 


-0.70 0.16 1.31 


j,D+-^K^l{1400)^^2^ 


0.58 -0.01 1.31 


j,D+^K° (1270) ^^2-^ 


7.12 -1.23 1.35 


_^D+^K°{1400)^^2^ 


-5.32 -0.87 1.35 



Table 3: Parameters appearing in the fit function for the form factors of the Dg 
Ki{1270)eu and Dq Ki{UOO)eu decays at = 8 GeV^, M| = 6 GeV^ and Ok, = 37°. 





a b trifit 




a b mfit 


^^-.^,(1270) ^^2) 


2.12 -0.30 1.27 


^^^i..(1400)^^2) 


-7.44 3.10 1.27 


^Z.^^.(m0)(^2) 


-1.52 1.10 1.37 


^i.-.K,(U00)(^2) 


2.70 -1.78 1.37 


^£)^ifi(1270)^^2^ 


0.27 -0.08 1.29 


y.D^i^l(1400)^^2j 


-0.75 0.26 1.29 


^D^i..(1270)(^2) 


1.68 -0.22 1.31 


^D^Ki(1400)|^^2^ 


-4.00 0.51 1.31 


^D+^i^0(1270)^^2-) 


1.29 0.93 1.30 


jD+^Kl{1400)^^2-^ 


-9.18 3.91 1.30 


^I3+-^i^0(1270)^^2^ 


-2.14 1.49 1.53 


^D+-.if0(1400)^^2) 


4.03 -2.59 1.53 


^I?+^K0(1270)^^2^ 


-0.45 0.14 1.32 


^D+^X0(1400)^^2^ 


0.79 -0.06 1.32 


^13+ -.Ki' (1270) ^^2^ 


4.78 -1.42 1.37 


^D+^X0(1400)^^2^ 


-7.57 -0.40 1.37 



Table 4: Parameters appearing in the fit function for the form factors of the Dq — > 
Xi(1270)£i/ and Dq Xi(1400)£i/ decays at Mf = 8 GeV^, M| = 6 GeV^ and 9k, = 58°. 
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a b rufit 




a b rufit 




5.48 -1.48 1.23 




2.81 -0.54 1.23 




-2.95 2.02 1.33 


y.D^i^i(1400)^^2j 


-1.71 1.22 1.33 




0.61 -0.17 1.25 


^£>-^Ki(1400)^^2^ 


0.35 -0.12 1.25 




3.90 -0.66 1.29 


^D^K.iim^^2^ 


2.10 -0.28 1.29 


j.D+-^K° (1270) ^^2^ 


7.23 -2.37 1.27 


^D+^i^0(1400)^^2) 


2.10 0.66 1.27 


j,D+^KO (1270) ^^2^ 


-4.27 2.83 1.48 


jD+^Kl{lA00)^^2-^ 


-2.43 1.67 1.48 


j,D+-.K0{1270) ^^2-^ 


-0.80 0.14 1.30 


^D+^if0(1400)^^2^) 


-0.54 0.15 1.30 


j,D+ (1270) ^^2^^ 


7.05 0.28 1.36 


^D+^X0(1400)^^2^ 


5.62 -1.44 1.36 



Table 5: Parameters appearing in the fit function for the form factors of the Dq 
Ki{1270)iiy and Dq Ki{UOO)eiy decays at Mf = 8 GeV\ M| = 6 GeV^ and Ok, = -37°. 

At the end of this section, we would hke to discuss the numeric values of the differential 
decay rates as well as the branching ratios for the considered semileptonic and nonleptonic 
transitions. 

semileptonic 

The dependence of the longitudinal and transverse components of the differential decay 
width for the semileptonic Dg Kiii/ decays is shown in Figs. 9-20 at 9ki = ±37°. In 
these figures, the total decay widths related to each decay are also depicted. To calculate 
the branching ratios of the semileptonic decays, we Integrate Eq. (30) over in the 
whole physical region and using the total mean life-time r^jo = 0.41 ps, Td+ = 1-04 ps 
and = 0.50 ps [23]. The values for the branching ratio of these decays are obtained 
as presented in Table 7. The errors in this Table are estimated by the variation of the 
Borel parameters and M|, the variation of the continuum thresholds Sq and s'q and 
uncertainties in the values of the other input parameters. 

nonleptonic 

For estimating the branching ratio of the nonleptonic D — >■ KiTT decay, first the values of 
the form factors at = ml are calculated as shown in Table 8. Inserting these values in 
Eq. (36) and using = 0.97377 ± 0.00027 [23], = 0.139 GeV and = 0.133 GeV, 
we obtain the values for the branching ratio of these decays as presented in Table 9. In 
comparison, we also include the experimental values and upper limits in this Table. This 
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a b rufit 




a b rufit 




3.86 -0.91 1.24 




4.88 -1.28 1.24 




-2.17 1.49 1.35 


y.D^i^i(1400)^^2j 


-2.57 1.80 1.35 




0.44 -0.10 1.26 


^£>-^Ki(1400)^^2^ 


0.56 -0.18 1.26 




2.97 -0.61 1.27 


^D^i.,(1400)^^2) 


3.38 -0.48 1.27 


j.D+-^K° (1270) ^^2^ 


5.73 -2.15 1.29 


^D+^i^0(1400)^^2) 


4.88 -0.48 1.29 


j,D+^KO (1270) ^^2^ 


-3.14 2.07 1.49 


jD+^Kl{lA00)^^2-^ 


-3.70 2.50 1.49 


j,D+-.K0{1270) ^^2-^ 


-0.58 0.08 1.32 


^D+^if0(1400)^^2^) 


-0.78 0.17 1.32 


j,D+ (1270) ^^2^^ 


4.69 0.71 1.35 


^D+^X0(1400)^^2^ 


7.84 -1.20 1.35 



Table 6: Parameters appearing in the fit function for the form factors of the Dq — > 
Ki{1270)iiy and Dq Ki{UOO)eiy decays at Mf = 8 GeV\ M| = 6 GeV^ and Ok, = -58°. 

Table shows that for the D° ii'r(1270)7r+, i^f (1400)7r+ and D+ K^{U00)7i+ 

cases, the different values of mixing angle give the values of branching ratios in good 
agreement with the experimental results but for — > 1270)77+ decay, the values of 
the branching ratios at different values of are about one order of magnitude more than 
that of the experimental expectation. 

In summary, we analyzed the semileptonic Dq — > Kiiu transition with q = u,d,s in 
the framework of the three-point QCD sum rules and the nonleptonic D — > Kiir decay 
within the factorization approach. We calculated Dq to iri(1270) and i^i(1400) transition 
form factors by separating the mixture of the ii'i(1270) and A'i(1400) states. Using the 
transition form factors of the D ^ Ki, we analyzed the nonleptonic D KiTT decay. 
We also evaluated the decay amplitude and decay width of these decays in terms of the 
transition form factors. The branching ratios of these decays were also calculated at different 
values of the mixing angle ■ For the non leptonic case, a comparison of the results for 
the branching ratios with the existing experimental results was also made. 
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01- 


37 58 -37 -58 


Br{D^ (1270)£i/) 


[3.59 ±0.29 1.03 ±0.10 5.34 ±0.21 2.84 ± 0.25] x lO'^ 


Br{D+ is:0(1270)£z/) 


[9.47 ±0.45 2.70 ±0.25 14.07 ± 1.22 7.57 ± 0.35] x lO'^ 


DiiDt Kl{m{})(p) 


[7.81 ±0.11 2.09 ±0.21 12.51 ±1.16 6.91 ± 0.32] x 10"^ 




[1.09 ±0.10 1.78 ±0.15 0.85 ±0.02 1.20 ± 0.11] x lO'^ 


Br{D+ K^{UOO)iu) 


[2.93 ±0.25 4.75 ±0.29 1.27 ±0.10 3.20 ± 0.27] x 10"=^ 


Br{D+ K^{Um)lv) 


[3.44 ±0.29 5.88 ±0.34 1.49 ±0.13 3.96 ± 0.29] x 10"^ 



Table 7: The values for the branching ratio of the semileptonic Dq — > i^i(1270)£i/ and 
Dq — > Ki{lAQ{))lv decays at different values of the Ok^. 





37 58 -37 -58 




37 58 -37 -58 


rD^Ki{l27Q) 

Jv 


3.24 1.82 4.04 2.95 


Jv 


-3.45 -4.42 2.30 3.65 


JO 


-0.73 -0.42 -0.91 -0.67 


rD^Ki{1400) 

Jo 


0.70 0.92 -0.47 -0.75 


nD-*Ki{1270) 
Jl 


0.34 0.20 0.45 0.32 


Jl 


-0.36 -0.49 0.25 0.41 


J2 


2.67 1.55 3.32 2.49 


h 


-2.81 -3.65 1.87 3.03 



Table 8: The values of the form factors of the D i^i(1270) and D Ki{UOO) for 
Mf — 8 GeV'^, M| — 6 GeV"^ at — ml and different values of the mixing angle 6ki- 





37 58 -37 -58 


Exp [23 


Br(D^ (1270)7r+) x lO'^ 


1.45 ±0.11 0.75 ±0.06 2.26 ±0.18 1.23 ±0.11 


1.15 ±0.32 


Br{D+ K^{127Q)tt+) x lO'^ 


3.75 ±0.29 1.23 ±0.10 5.85 ± 0.37 3.18 ±0.25 


< 0.7 


Br(D^ irr(1400)7r+) x 10"^ 


0.60 ±0.04 1.00 ±0.12 0.26 ±0.02 0.73 ±0.04 


< 1.2 


Br{D+ is:0(1400)7r+) x lO'^ 


2.57 ±0.21 3.63 ±0.31 1.71 ± 0.13 2.78 ± 0.24 


3.8 ± 1.3 



Table 9: The branching ratios of the nonleptonic D iri(1270)7r and D iri(1400)7r 
decays at different values of 9ki- 
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Figure 5: The dependence of the form factors on 9ki at = for D —>■ Ki{ynQi)£i/ decay. 
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Figure 6: The dependence of the form factors on Ok^ at = for D — > Ki{lAQ{))lv decay. 
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Figure 7: The dependence of the form factors on Ok^ at = for Dg —>■ Ki{ynQi)£v decay. 
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Figure 8: The dependence of the form factors on Ok^ at = for Dg — > Ki{lAQ{))lv decay. 
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Figure 9: The dependence of the dVtot/dq^, dVT/dq^ and dVL/dq^ on at Ok-^ — 37° for 
(1270)^1/. 




Figure 10: The dependence of the dTtot/dq^, dVT/d(f and dV^/dq^ on q"^ at = 37° for 
D+ i^0(1270)£z/. 
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Figure 11: The dependence of the dTtot/dq^, dVT/dq^ and dVL/dq^ on at Ok^ — 37° for 




Figure 12: The dependence of the dTtot/dq^, dVT/d(f and dV^/dq^ on at = 37° for 
^ AV(1400)£t/. 
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Figure 13: The dependence of the dTtot/dq^, dVT/dq^ and dVL/dq^ on at Ok^ — 37° for 
^ 1400)^1/. 




Figure 14: The dependence of the dTtot/dq^, dVT/d(f and dV^/dq^ on at = 37° for 
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Figure 15: The dependence of the dTtot/dq^, dVT/dq^ and dVL/dq^ on at Ok^ — —37° for 




Figure 16: The dependence of the dVtot/dq^-i dTx/dq^ and dT^/dq^ on q"^ at 9ki = —37° for 
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Figure 17: The dependence of the dTtot/dq^, dVT/dq^ and dVL/dq^ on at Ok^ — —37° for 




Figure 18: The dependence of the dVtot/dq^-i dTx/dq^ and dT^/dq^ on q"^ at 6'^^ = —37° for 
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Figure 19: The dependence of the dTtot/dq^, dVT/dq^ and dVL/dq^ on at Ok^ — —37° for 
D+ Kl{imQ)£v. 




Figure 20: The dependence of the dVtot/dq^-i dTx/dq^ and dT^/dq^ on q"^ at 9ki = —37° for 
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